Emerging infectious diseases and animal social systems

Emerging infectious diseases threaten a wide diversity of animals, and important questions remain concerning disease emergence in socially structured populations. We developed a spatially explicit simulation model to investigate whether—and under what conditions—disease-related mortality can impact rates of pathogen spread in populations of polygynous groups. Specifically, we investigated whether pathogen-mediated dispersal (PMD) can occur when females disperse after the resident male dies from disease, thus carrying infections to new groups. We also examined the effects of incubation period and virulence, host mortality and rates of background dispersal, and we used the model to investigate the spread of the virus responsible for Ebola hemorrhagic fever, which currently is devastating African ape populations. Output was analyzed using regression trees, which enable exploration of hierarchical and non-linear relationships. Analyses revealed that the incidence of disease in single-male (polygynous) groups was significantly greater for those groups containing an average of more than six females, while the total number of infected hosts in the population was most sensitive to the number of females per group. Thus, as expected, PMD occurs in polygynous groups and its effects increase as harem size (the number of females) increases. Simulation output further indicated that population-level effects of Ebola are likely to differ among multi-male–multi-female chimpanzees and polygynous gorillas, with larger overall numbers of chimpanzees infected, but more gorilla groups becoming infected due to increased dispersal when the resident male dies. Collectively, our results highlight the importance of social system on the spread of disease in wild mammals.     

Molecular dynamics simulation of isothermal crystallisation of polymer chains around single polymer lamella

The isothermal crystallisation of polyethylene (PE) chains around single PE lamella in vacuum is investigated by molecular dynamic simulation. The crystallisation process is analysed in terms of the orientational order parameters, principal moments of inertia for the simulated systems. The effects of charge interactions between the polymer chains and lamella are discussed. It is found that the crystallisation process for uncharged systems can be divided into three stages: (1) adsorption, (2) orientation and (3) arrangement. The single polymer lamella changes a little during the three stages. PE chains are arranged parallel to the chain direction of the stems in the crystalline state. When considering the effect of charge interactions between the polymer chains and lamella, a different crystallisation process appears. The single polymer lamella is affected by the charged polymer chains.     

How good are comparative models in the understanding of protein dynamics?

The 3D structure of a protein is essential to understand protein dynamics. If experimentally determined structure is unavailable, comparative models could be used to infer dynamics. However, the effectiveness of comparative models, compared to experimental structures, in inferring dynamics is not clear. To address this, we compared dynamics features of ~800 comparative models with their crystal structures using normal mode analysis. Average similarity in magnitude, direction, and correlation of residue motions is >0.8 (where value 1 is identical) indicating that the dynamics of models and crystal structures are highly similar. Accuracy of 3D structure and dynamics is significantly higher for models built on multiple and/or high sequence identity templates (>40%). Three-dimensional (3D) structure and residue fluctuations of models are closer to that of crystal structures than to templates (TM score 0.9 vs 0.7 and square inner product 0.92 vs 0.88). Furthermore, long-range molecular dynamics simulations on comparative models of RNase 1 and Angiogenin showed significant differences in the conformational sampling of conserved active-site residues that characterize differences in their activity levels. Similar analyses on two EGFR kinase variant models highlight the effect of mutations on the functional state-specific αC helix motions and these results corroborate with the previous experimental observations. Thus, our study adds confidence to the use of comparative models in understanding protein dynamics.     

A Refined Open State of the Glycine Receptor Obtained via Molecular Dynamics Simulations

1. Download : Download high-res image (389KB)
2. Download : Download full-size image

     

Identification and in silico prediction of metabolites of tebufenozide derivatives by major human cytochrome P450 isoforms

Cytochrome P450 (CYP) enzymes constitute a superfamily of heme-containing monooxygenases. CYPs are involved in the metabolism of many chemicals such as drugs and agrochemicals. Therefore, examining the metabolic reactions by each CYP isoform is important to elucidate their substrate recognition mechanisms. The clarification of these mechanisms may be useful not only for the development of new drugs and agrochemicals, but also for risk assessment of chemicals. In our previous study, we identified the metabolites of tebufenozide, an insect growth regulator, formed by two human CYP isoforms: CYP3A4 and CYP2C19. The accessibility of each site of tebufenozide to the reaction center of CYP enzymes and the susceptibility of each hydrogen atom for metabolism by CYP enzymes were evaluated by a docking simulation and hydrogen atom abstraction energy estimation at the density functional theory level, respectively. In this study, the same in silico prediction method was applied to the metabolites of tebufenozide derivatives by major human CYPs (CYP1A2, 2C9, 2C19, 2D6, and 3A4). In addition, the production rate of the metabolites by CYP3A4 was quantitively analyzed by frequency based on docking simulation and hydrogen atom abstraction energy using the classical QSAR approach. Then, the obtained QSAR model was applied to predict the sites of metabolism and the metabolite production order by each CYP isoform.     

Hard-sphere/pseudo-particle modelling (HS-PPM) for efficient and scalable molecular simulation of dilute gaseous flow and transport

Continuum methods are not accurate enough for flows at high Knudsen numbers, whereas rigorous molecular dynamics (MD) methods are too costly for simulations at practical dimensions. Hard-sphere (HS) model is a simplified MD method efficient for dilute gaseous flow but is of poor parallelism due to its event-driven nature, which sets a strong limitation to its large-scale applications. In this work, pseudo-particle modelling, a time-driven modelling approach is coupled with HS model to construct a scalable parallel method capable of simulating flows and transport processes at high Knudsen numbers without losing necessary molecular details in describing their macro-scale behaviours. The method is validated in several classical simulation cases and its performance is evaluated to be favourable. To demonstrate the potential applications of this method, we also simulate the diffusion of small molecules in multi-scale porous media which is related to catalysis, material preparation and micro chemical engineering in the long term.     

Atomic insight into prion disorder: An intricate detail gained by 0.5 μs molecular dynamics simulation of preventive G127V and deleterious D178V mutation in prion protein

In this study we are looking into two contradicting mutations found in prion protein (PrP) viz G127V and D178V, that are reportedly protective and pathogenic, respectively. Despite significant advances in comprehension of the role of pathogenic mutations, the role of protective mutation in amyloid fold inhibition still lacks a substantial basis. To understand the structural basis of protective mutation, molecular dynamics simulation coupled with protein-protein docking and molecular mechanics/Poisson-Boltzmann surface area analysis was used to understand the instant structural variability brought about by these mutations alone and in combination on PrP and prion-prion complex. Atomic-scale investigations successfully revealed that the binding pattern of prion-prion varies differentially in protective and pathogenic mutations with secondary structure showing distinct contrasting patterns, which could supposedly be a critical factor for differential prion behavior in protective and pathogenic mutations. Considering the reported role of an amyloid fold in prion-prion binding, the contrasting pattern has given us a lead in comprehending the role of these mutations and has been used in this study to look for small molecules that can inhibit amyloid fold for prion-prion interaction in pathogenic mutant carrying PrP.     

A sequential test for assessing observed agreement between raters

Assessing the agreement between two or more raters is an important topic in medical practice. Existing techniques, which deal with categorical data, are based on contingency tables. This is often an obstacle in practice as we have to wait for a long time to collect the appropriate sample size of subjects to construct the contingency table. In this paper, we introduce a nonparametric sequential test for assessing agreement, which can be applied as data accrues, does not require a contingency table, facilitating a rapid assessment of the agreement. The proposed test is based on the cumulative sum of the number of disagreements between the two raters and a suitable statistic representing the waiting time until the cumulative sum exceeds a predefined threshold. We treat the cases of testing two raters' agreement with respect to one or more characteristics and using two or more classification categories, the case where the two raters extremely disagree, and finally the case of testing more than two raters' agreement. The numerical investigation shows that the proposed test has excellent performance. Compared to the existing methods, the proposed method appears to require significantly smaller sample size with equivalent power. Moreover, the proposed method is easily generalizable and brings the problem of assessing the agreement between two or more raters and one or more characteristics under a unified framework, thus providing an easy to use tool to medical practitioners.